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Abstract 

We study the shear mode in the gauge/gravity correspondence at 
finite temperature. First, we confirm the general formula for the shear 
viscosity in an arbitrary background metric which includes a black 
hole in the fifth dimension. We then derive a general formula for the 
shear mode relaxation time which appears in the theory of relativistic 
dissipative fluid dynamics; it agrees with known expressions in the 
limit of conformal fields. These results may be useful in relativistic 
viscous fluid descriptions of high energy nuclear collisions at RHIC 
and LHC. 
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1 Introduction 

The AdS/CFT correspondence p] or, more generally, the relationship be- 
tween gauge theories and higher dimensional gravity theories, is a useful tool 
for studying strongly coupled plasmas, such as the fluid created at RHIC 
(Relativistic Heavy Ion Collider) during heavy ion collisions. One can extract 
hydrodynamic transport coefficients of the plasma by examining perturba- 
tions in the dual gravity theory. Usually this is done by either invoking the 
AdS / CFT prescription to compute correlation functions of the stress-energy 
tensor [21 El HI El E] or by applying appropriate boundary conditions and 
examining the resulting dispersion relation 0E]. In [H], the shear viscosity 
of currents defined on a black hole's stretched horizon was computed for a 
general class of supergravity backgrounds. The value of the stretched hori- 
zon shear viscosity agrees in all cases with the shear viscosity computed via 
the methods previously listed [10]. Though this agreement is not understood 
at present, this membrane paradigm approach provides yet another method 
for computing the transport coefficients of the strongly coupled plasma (as- 
suming that the agreement between the membrane paradigm and AdS / CFT 
approaches is absolute). We employ the membrane paradigm approach in 
this work. 

Hydrodynamics is an effective theory which describes the dynamics of 
a thermal system on length and time scales much larger than any micro- 
scopic scale. Here we consider perturbations of the fluid in the hydrody- 
namic regime; in this case, the energy and momentum of the perturbations 
must be much less than the temperature of the plasma u,q <C T. In a hy- 
drodynamic theory, the spatial components T iJ of the stress-energy tensor 
are constructed from the conserved quantities T 00 = e (energy density) and 
rpOi = n t ( m omentum density). These 'constitutive relations' arise by writing 
all possible terms consistent with the symmetries present, taking into account 
equilibrium thermodynamics, and by ignoring higher derivative terms since 
the hydrodynamic modes are assumed to be slowly varying in space and time. 

As an example, consider a translationally invariant, four dimensional 
theory in flat space. To first order in the perturbations and derivatives the 
constitutive relations are 

T ij = 5 ijrp + v 2 5 A _ l^TT* _1( + gj^i _ ls ij d k 7T k ) . (1) 

w w V 3 J 

Here the Latin indices i, j run over the 3 spatial coordinates. The symbols P 
and w represent the equilibrium pressure and enthalpy density, respectively, 
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and v s , r\ and ( represent the speed of sound, shear and bulk viscosities. Of 
course the energy-momentum tensor obeys the conservation equation, with 
the Greek indices fi, v running over all four coordinates, 

dfjT^ = . (2) 

These linearized hydrodynamic equations admit two normal modes, cor- 
responding to whether the momentum density fluctuations are transverse or 
longitudinal to the fluid flow. Transverse fluctuations lead to the shear mode 
which is the one considered in this paper. There is also the sound mode (from 
longitudinal momentum fluctuations), and the diffusive mode (in the pres- 
ence of a conserved current); these modes are not considered in this work. 
For a more complete introduction to relativistic hydrodynamics see [TTj . 

Hydrodynamic fluctuations in the gauge theory plasma are dual to per- 
turbations of a specified gravitational background in the dual theory. Through- 
out this paper we use an effective 5-dimensional metric. In principle, one can 
view this metric as being derived by dimensionally reducing some theory in 
a higher dimension. Our results are easily generalizable to different numbers 
of dimensions in the same way as the results of [9]. Our convention is to 
use t, x, y, z to denote the usual four space-time coordinates; the coordinate 
r denotes the extra dimension. 

In the case of the shear mode, one only needs to consider metric pertur- 
bations g^ u — > g^ u + h^ u . The components of the perturbation h^ u can be 
decomposed into three irreducible sets; the field equations for each set de- 
couple from the others. The sets are found by classifying the perturbations 
under 0(2) rotations [I]; see also [8]. The set of gravitational perturbations 
dual to the hydrodynamic shear mode have h y Q, h yz , and h yr all nonzero with 
all other components of h^ u vanishing. We assume a standard gauge choice 
h^ r — which leaves us with only two nonvanishing components of h^ u . 

The hydrodynamic transport coefficients are found by solving the lin- 
earized Einstein equations for the perturbations with appropriate boundary 
conditions. The resulting dispersion relation for these gravitational perturba- 
tions can be compared with the expected dispersion relation from the bound- 
ary (gauge) hydrodynamics. In first order hydrodynamics, this method has 
been used to calculate the speed of sound and the shear and bulk viscosities 
in a wide variety of gravity duals. In [9J, the authors also derive a gen- 
eral formula for the shear viscosity which depends only on components of 
the dual metric; this formula is applicable for a large class of gravitational 
backgrounds. 
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Recently, several groups have begun to extend these computations to 
second order hydrodynamics [121 H31 EH]- Currently, there is no unique for- 
malism for second order hydrodynamics, though a frequently used one is the 
"Miiller-Israel-Stewart" formalism [15J. All such extensions of first order hy- 
drodynamics attempt to repair problems that the first order theory has with 
causality, and necessarily introduce a set of new transport coefficients. One 
such quantity is the "shear relaxation time" r, and appears in the next to 
leading order in the shear dispersion relation, as explained in Sec. 5 below. 

In this paper we derive a general formula for the shear relaxation time 
analogous to the formula for the shear viscosity given in [9j. This formula is 
applicable to a wide variety of gravitational duals, and reproduces all known 
results in the literature of r. 

This paper is organized as follows. In Sec. 2 we define the relevant 
gravitational background, the background Einstein equations, and the con- 
servation equations that must be satisfied by the perturbed stress-energy 
tensor. Our general relativistic conventions are those of [16]. In Sec. 3 we 
set up the gravitational shear perturbations, and determine the necessary 
linearized Einstein equations and boundary conditions that must be satisfied 
by the perturbations. In Sec. 4 we proceed to solve these equations to lowest 
nontrivial order in q. In doing so, we reproduce the formula for the shear 
viscosity which was derived by a different method in [9]. In Sec. 5 we extend 
the calculation to the next order in q, thus deriving the corresponding general 
formula for the relaxation time r. In Sec. 6 we consider applications of this 
result to certain phenomenological models of QCD, namely, the "Hard Wall" 
and "Soft Wall". We make some conclusions and mention some prospects 
for further work in Sec. 7. 

2 Metric Perturbations 

We consider a background metric of the generic form 

ds 2 = g 00 (r)dt 2 + g xx (r)\dx\ 2 + g„{r)dr 2 (3) 

where the three functions goo-, 9xx and g rr depend on the fifth dimensional 
variable r only. Our convention is such that goo = —1 and g xx = 1 in 
flat 4-dimensional space. The resulting Riemann-Christoffel curvature tensor 
and Ricci tensor are given in the appendix. We assume that there is a 
horizon located at r = r . As the horizon is approached we assume that the 
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metric behaves as follows: g„ — > r y r /{r — r ) and g o — * ~lo( r ~ r o) while 
g xx approaches the finite value g xx (r )- Thus y/^g — > \/loJr9xx( r o) ■ The 
Hawking temperature is 

^ (4) 

The matter and radiation fields which generate this metric are not known. 
Nevertheless, the energy-momentum tensor which generates it can be inferred 
from the Einstein equations. They can be written either as 

Guv = Rfiu — = —^TrGdT^ (5) 

or as 

R, u = -8nG d (r^ - ^ T W) (6) 

where the number of space-time dimensions is d and Gd is the corresponding 
Newton constant. In this paper we focus on d — 5. With the notation used 
in the appendix this yields explicit expressions for the background energy- 
momentum tensor, 

Too = : ttt^ (—To + 3F X + F r ) , 
16nG 5 

T xx Tyy T zz „ (To -\- F x -\- T r ) , 

107TG5 

T rr = TTTT (To + 3F X - F r ) . (7) 

107T(jr5 

For a space with constant negative curvature Fq = F x = F r = 4/L 2 and 
R = 20/L 2 . Then T„ v = (3/47rG 5J L 2 )<7 MI ,. 

For metric perturbations g^ u — > g^u( r ) + h^ u (t,x.,r) we allow for a per- 
turbation in the energy-momentum tensor which is first order in and its 
derivatives. 

= Tffi + Tffi. (8) 

Here satisfies the Einstein equations with the background metric g^ v 
given in Eq. (j3J), see Eq. (J7J), and T?X is its correction to first order in 
hfa,. In addition it must satisfy the equation for conservation of energy and 
momentum. Expanding the exact equation 

T>" = (9) 
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to first order one gets 

Tfi* + T (o) + m; x = o (10) 

where (<5r)^ A is the change in the affine connection to first order in h^ v . For 
a space of constant negative curvature, as in a conformal field theory with a 
gravity dual, the answer is easy to obtain. It is just Tj^ = (3/ 'AttG^L 2 )^ 1 ' . 
In the general case it is more involved. 



3 Shear Modes 



We consider a metric fluctuation of the following form: 



h y0 (z,t,r) 
hy Z (z,t,r) 



g xx (r)A(r)e 
9xx(r)B(r)e 



i(qz-u)t) 
i(qz—U)t) 



;i2) 



Fluctuations in the metric are intimately connected with fluctuations in the 
energy-momentum tensor density. For linear response it is natural to assume 
that Tj^ has the same harmonic variation in z and t. The conservation 
equations then take the following form: 



1 



d 



99oo dr 
1 d 



o, 



-99xx dr 



1 



d 



-99xx dr 



Oz 

(i) 



-iu;T$g™h 0y + iqT$g™h zy 
1 d 



99ooT {l) j 
gg xx Tfi)=0, 

99xxT(f^J 
0, 



99xx dr 



(V^ggxxT^) = 0, 



1 d / \ 1 dg 00 00 

tvH7<7rrT (1) j - 5 </ — T (1) 



1 rr ^7xx 



(13) 
(14) 



(15) 
(16) 



(7^ + 7$ + 7$)=0. (17) 
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These equations are naturally satisfied by the choice 



T (i) 



16ttGf 



(Fq + F x + F r ) h 



fJbV 



(18) 



for arbitrary values of u and q. This is a particular solution of the inhomoge- 
neous equations, that is, for given h y0 and h yz . To it one may add solutions 
to the homogeneous equations, but this is not necessary to obtain the shear 
dispersion relation. 

The above analysis leads to three nontrivial Einstein equations of motion, 



d(1) 
^yO 

BP 

L yr 



F x h y o , 
F h 
0. 



(19) 
(20) 
(21) 



They can be expressed in terms of the unknown functions Air) and B(r), 



1 d 
'—g dr 

1 



— rr 00 

-99 9 9xx 



dA 
dr 



d 



-g dr 



-99 



ug 



00 



, r dB_ 
dr 
dA 



= qg m {qA + uoB) 

ug 00 (qA + loB) , 
„„dB 



qg 



(22) 



(23) 



(24) 



dr dr 

Of course these cannot all be independent since there are only two functions 
A and B. Substitution of Eq. (01]) into Eq. ([23]) reproduces Eq. for 
example. 

A boundary condition is needed near the horizon. From the structure of 
the equations near the horizon we observe that A ~ (r — r )B as r — > r . In 
this limit Eq. fl23]) implies that 



70/ sd 

-y- ro) dr- 



(25) 



which has two independent solutions corresponding to waves incoming to the 
horizon and outgoing from the horizon, 



h 



fin [t + J — ln(r - r ) + /out * - W — ln(r - r 



7o 



(26) 



The boundary condition is to select only the incoming waves. This results in 
the following requirement: 



.&r-r„)^ as r^r . (27) 
at y 7 r ar 

For the exponential time dependence assumed in Eq. (jT2l) this is 

- iuhy Z = J — (r - ro)—^- as r -> r . (28) 



The idea proposed in [9] is to apply this condition at a distance such that 
r h~ r o ?"o, which is referred to as a stretched horizon. This is then inserted 
on the left side of Eq. (|23|) to obtain 

(r-r )^- = -iJ—{qA + uB) as r->r h . (29) 
dr y 70 

The reason for introducing the stretched horizon has to do with the dispersion 
relation uj = —iDq 2 + ■ ■ -. Although \A\ \B\ as r — > r , consideration of 
very small but nonvanishing wave- number q will result in q\A\ ~ \u)B\ and 
eventually q\A\ 3> |^-B| for a fixed value of r as g — > 0. The latter can be 
satisfied self-consistently if — r$ is not too small. A more precise condition 
will be given below. 



4 Shear Viscosity 

The shear modes are transverse and strongly over-damped. At low frequency 
and wave-number they have the dispersion relation to = —iDq 2 where D is 
a diffusion constant. In terms of the shear viscosity r\ and enthalpy density 
w = Ts = P + p it is D = rj/w. Upon inspection, and with the expectation 
of a diffusive mode, we expand the functions in powers of q as follows: 

oo 

A(r) = E^ (2n) (0, (30) 

n=0 

oo 

B(r) = f^q 2n+1 B^ n+1 \r). 

n=0 
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From Eqs. (1221) and (1231) we find the lowest order solutions 



#(r) = a + a 1 f dr' ^9^) ^ 
Jo ° V -g(r')g xx {r>) 

B U(r) = bp + h f dr' 9 ; r{rl) . (32) 
J oo . I — g(r') 



Substitution into Eq. (JMD, together with u = —iDq 2 , gives the relationship 
D = %b\ja\. The Dirichlet boundary condition on A and B at infinity gives 
a = b = 0. 

It is now a straightforward matter to substitute the lowest order solution 
into the boundary condition on the stretched horizon to obtain the shear 
diffusion constant, 

D = - ir 7*'^' . (33) 

V-^oo(^o)fl'rr('"o) Jr ° y-g(r)g xx (r) 

Here the limit — > ro has been taken without any harm. This is identical 
to the result obtained in [9j in which Maxwell's equations were solved (in the 
specified metric) rather than Einstein's equations. 

A frequently used metric is one that is dual to M = 4 SU(iV) supersym- 
metric Yang Mills theory at finite temperature and in the limit N — > oo and 
g 2 N -> oo, 



ds2 = v 



2 _ L 2 



'" f(r)dt 2 + \dx\ 2 +-^—dr 2 . (34) 
j r 2 f(r) 



Here f(r) = 1 — (ro/r) 4 . Using this metric yields 

D 



AttT 
i] 1 

s 4tt ' 



(35) 



which are well-known results. 

The self-consistency of the assumption that \uB\ could be neglected com- 
pared to q\A\ near the stretched horizon can now be checked. It requires the 
following: 

( 4ttT\ . , 

r exp I —jjz J < r h - r < r . (36) 
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Obviously there is a finite range of q for which this condition can be met. 
Nevertheless the derived expression for the diffusion constant and the ratio 
of shear viscosity to entropy density remains valid. 



5 Shear Relaxation Time 

The lowest order dispersion relation for shear modes is u — —iDq 2 . The 
partial differential equation which describes such a mode is the diffusion 
equation 

§£*- D &* = (37) 

where tp represents a component of the shear-stress tensor. The frequency 
is obviously an expansion in even powers of the wavenumber. The next 
correction to the shear dispersion relation given above is fourth order in q. 
If we want to write a partial differential equation to describe the shear mode 
we should add a term proportional to the fourth spatial derivative of the 
field. For consistency, there ought to be a term proportional to the second 
temporal derivative of the field as well. Therefore we are led to consider the 
partial differential equation 

t 4^^'- d ^ + ^ d2 &^° (38) 

where two time constants T\ and r 2 have been introduced. This leads to a 
shear dispersion relation 



iDq 2 (l + rDq 2 + • • •) (39) 



where only the combination r = T\ + r 2 appears. We refer to r as a "shear 
relaxation time" because it is associated with the shear mode and it has the 
dimensions of time. The term is somewhat of a misnomer because it does not 
imply that fields falloff as exp(— t/r); see Sec. 6. In [12] the quantity which 
we denote as r was denoted rn and was believed to coincide with the "re- 
laxation time" in the Muller-Israel- Stewart formulation of second order fluid 
dynamics [15] . However, in [13] it was argued that the quantity we denote as 
r contains not only the Miiller-Israel-Stewart rn, but also contributions from 
(currently unformulated) third order fluid dynamics. Actually this is quite 
apparent from Eq. fl38|) . In this way the authors of [13] explained the ap- 
parent discrepancy that arises when inferring rn from the shear mode versus 
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the sound mode. To clarify, we emphasize that our r is the same quantity 
that is computed in [T2], but is not the same as the r n of Miiller, Israel and 
Stewart. However, it is necessary input for the higher order relativistic dis- 
sipative fluid dynamic description of the flow of matter and radiation. We 
shall now derive a formula for r. 

The lowest order solution, which was all that was needed to determine 
the shear viscosity, was 



dr 



/ -goo(r')g rr (r') 



^-9(r')9xx(r') 



poo 

B w (r) = -iDax / dr 



, 9rr(r') 



(40) 



Substituting these into Eqs. (J22T) - f|23T) leads to the next higher order terms, 
, grr(r')g 00 (r') 



A (2) (r) = a 2 + dr 



\l-g{r')gxx{r') 



B^\r) = b 2 + [ r dr 

J oo 



-9(f) 



-\a 3 + [ dr"J-g(r")g 00 (r")A^(r") 

l\ I Joo 

iD f dr"^-g(r")g°°(r")A^(r")) 



(41) 

Dirichlet boundary conditions at infinity require a 2 = b 2 = 0. Finally, sub- 
stitution of these solutions into Eq. (1241) leads to 



D 2 r 



-D h % — . 

a\ ai 



(42) 



To complete the calculation the boundary condition on the stretched 
horizon (1291) must be satisfied. To this order it is 



dB^ 



(r - r J 



A^-iDB^) as 



dr 4irT 
This precisely determines the shear relaxation time 



(43) 



T 



-g(r ) 



\/-goo(r )g r r(r ) Jr ° yf-g(r) 



°° dr Jrr(r) 



1 



' Dirf 



(44) 
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where 

\/~9{t) [°° , , -gm{r')g rr {r') 



B{r) = v = , / dr' auuv J » rrv ' . (45) 
y-9oo(r)g r r(r) Jr y -g{r')g xx {r') 

In deriving Eq. (|44|) an integration by parts was performed to reduce the 
number of nested integrals. Note that D(r ) is just the diffusion constant D. 
The limit — > r has been taken since the integrals converge. 

The formula given in Eq. (j4"4l) is the main result of this paper. For the 
special metric (1541) it gives 

1 - In 2 

T = "ST" (46) 

which was already derived in p^2l [13] . In the literature this metric is referred 
to as AdS5(D3) for 5-dimensional anti-deSitter D3-brane. The shear relax- 
ation time has been computed for other conformal theories as well, including 
AdS4(M2) and AdS7(M5). The formula (1141) reproduces all of these known 
results. 



6 Hard and Soft Wall Models of QCD 

There are several phenomenological models for QCD that incorporate con- 
finement. The simplest of these is the hard wall model. A better one is the 
soft wall model as it represents the linear radial Regge trajectories for the 
families of p and a\ mesons quite well. Here we compute the ratio of shear 
viscosity to entropy density and the shear relaxation time for these models 
since they may be useful in applications to high energy nuclear collisions at 
RHIC and LHC. 

6.1 Hard wall model 

The hard wall model [17J uses the metric of Eq. (13~4|) but with r-space 
cutoff at some value r min . This leads to radial excitations of the vector and 
axial- vector meson spectra with the mass being linear in the radial quantum 
number n. Introduction of the new scale r min implies that the transport 
coefficients do not scale dimensionally with temperature T alone. 

There is a significant difference depending on whether r min < r or r min > 
tq. In the former case the introduction of a hard cutoff does not influence 
the calculations of the previous sections since it lies inside the black hole. 
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In the latter case it would seem to be indeterminate because the horizon 
lies beyond the wall and the boundary condition near the horizon cannot 
be implemented. On the other hand, we can stretch the boundary out from 



tq to r 



if the distance rv 



— To is so small that we can use the 
approximations goo( r ) ~ ~7o( r — r o) and g rr (r) ~ j r /(r — r ). Then the 
boundary condition on the stretched horizon which allows only for incoming, 
and not outgoing, waves can be applied. The answers are the same as Eqs. 
( 133]) and fpf4|) with r replaced by r min . One finds that 

1 

Arc 



1 = — for T > T r 



(47) 



and 



for T <T r 



(48) 



i] _ 1 T 

where the Hawking temperature is T = tq/tiL 2 and we have defined a critical 
hardwall temperature T c = r min /7rL 2 . If we are willing to tolerate a 5% de- 
viation of the metric from the limiting form given above then the expression 
(|48]) is only a good approximation for 0.9 < T/T c < 1. For lower tempera- 
tures it is difficult to cleanly separate the incoming from the outgoing waves. 
The hard wall model is too crude to access lower temperatures. Nevertheless 
it does give an indication of how rj/s might deviate from its conformal value 
when a confinement scale is introduced. Note that it does fall below the 
conjectured lower bound for rj/s of 1/47T when T < T c . 
The hard wall result for the shear relaxation time is 

1 - In 2 



T 



2nT 



for T>T r 



(49) 



and 



AttT 



^ 2 



1+ ^ 









In 








1- rr 



(50) 



We expect this result to be reasonably accurate only when 0.9 < T/T c < 1 
for the reasons mentioned above. In that case the formula can be expanded 
just below T c as follows: 



2nT 



In 2 



T r -T 



In 



8T C 
T c -T 



(51) 
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Note that the effect of the confinement scale is to reduce the shear relaxation 
time compared to the conformal limit. Numerically, the value of T c can be 
estimated as follows. The mass of the p-meson in the hard wall model is found 
by solving a particular wave equation whose solutions are Bessel functions. A 
boundary condition requires that Jo(m p L 2 /r min ) = from which one deduces 
that T c « rrip/7.556 ~ 102 MeV. This is just a characteristic temperature, 
not the critical temperature of a phase transition, since one may define it 
somewhat differently |18j . 



6.2 Soft wall model 

The soft wall model [T5] was developed to improve upon the hard wall model. 
In particular, it leads to linear Regge trajectories wherein the radial excita- 
tions of the vector and axial-vector meson spectra have mass-squared being 
linear in the radial quantum number n, in substantial agreement with data. 
It can be obtained by adding a dilaton field to the usual AdS metric. The 
piece of the action that is relevant for the meson mass spectra is 

Smeson = ~T~ / ^V^T^ + F|) (52) 

where <fi is the dilaton field and F£ u and Fj^ are the field strength tensors for 
the left and right handed gauge fields. In our coordinates, the dilaton profile 
which leads to the Regge behavior is <p(r) = cL 4 /r 2 , where c is a constant 
which can be determined by fitting the meson spectrum. If we compute 
the quantities 77/s and r, using the soft wall model of [T5], the results are 
the same as Eqs. (1551) and (14"6T) because the metric is exactly AdS 5 . There 
is an alternative formulation of the soft wall. Instead of adding a nontrivial 
dilaton, one keeps the dilaton constant while deforming the metric away from 
AdS§. This version of the soft wall model has been studied in 
The deformed metric is 



ds 2 = e~ h ^ 



y»2 



L2 (-f{r)dt 2 + \d*\ 2 ) + —dr 2 



(53) 



In order to preserve the essential feature of the soft wall model, namely, the 
linear Regge trajectories, one must choose 

h(r) = 20(r) (54) 
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so that the action (1521) is the same with regard to the mesons. We stress 
that this is a different implementation of the soft wall model than that of 
[T9] ; it is not simply a transformation from the string frame to the Einstein 
frame, because the dilaton profile is modified as well. Both pictures lead to 
linear Regge trajectories, but other physical quantities may differ, as we now 
demonstrate by computing rj/s and r. 
Using this metric in Eq. f|55|) yields 



2ttx(] 



1 + — (e- X0 - 1 
x v 



(55) 



where, following [IE] (see also 
via T c 2 = c/4 and where 



x = 12 



), we define a characteristic temperature 



In the high temperature limit 



rj 1 

S 47T 



while in the low temperature limit 



1 
s 



TV 

24 



TFT 



2 fT s4 



12 VtJ + ' 



(56) 



(57) 



(58) 



Similar to the hard wall model, the conformal limit of 1/47T is approached 
from below. See Fig. 1. Unlike the hard wall model, the dependence of 77/s 
is smooth and well-defined for all temperatures. 

One may be surprised that the shear viscosity is not greater than or equal 
to 1/47T given the proof of [24] . The reason is that [24] assumes a background 
metric with the property that F (r) = F x (r). Instead, in this version of the 
soft wall model, one finds by explicit computation that 



F (r) 



-3h'[ 



r e 



h(r) r 



r 3 L 2 



(59) 



which is nonzero as long as h(r) is not constant. Thus, we would only expect 
the general proof of [24J to apply in the case where h(r) is constant where the 
metric is pure AdS^. The fact that Fo(r) 7^ F x (r) which implies rj/s 7^ 1/4-7T 
may be a sign that a metric such as ( 1531) cannot be generated by conventional 
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supergravity matter fields. However, as is customary, we view this metric as a 
phenomenologically motivated effective 5D theory, and thus will not concern 
ourselves with its origin in this paper. 

The relaxation time can also be computed in the soft wall model. Ap- 
plying the formula (1440 we find the relaxation time to be 



x xe * 
ax 



2ttT Jo Xq — x 2 



1 




x + e x — 1 
x + e~ x ° - 1 



(60) 



Although this integral may be expressed in terms of special functions, the 
result is not very enlightening. The relaxation time is plotted in Fig. 2. As 
for the hard wall model, the conformal limit is approached from below. It 
is noteworthy that the shear relaxation time is negative for T/T c < 0.686. 
Whether this is unphysical or not is discussed in the next subsection. 
A high temperature expansion is easily derived. 
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(61) 



This expansion is accurate to better than 10% when T/T c > 0.9. 

Numerically, the value of T c can be estimated as follows. The mass of the 
p-meson in the soft wall model is found by solving a particular wave equation 
whose solutions involve Laguerre polynomials. A boundary condition leads 
to the relation m p = l^fc. The Hawking-Page analysis of the phase transition 
in [18] gives T c = m p /A w 192 MeV. A similar analysis in the hard wall model 
yields a transition temperature of approximately m p /6.354 [18] . 

6.3 Sign of the shear relaxation time 

Is a negative shear relaxation time unphysical? First consider real, positive 
q and complex u. A shear plane wave propagating in the ^-direction has the 
form 

t/j = ^p e^*-"*) = e ^-<^ V 7 ' (62) 

where the frequency has been decomposed into its real and imaginary parts 
as to = + ioJi- The real part vanishes and the imaginary part is 

uj = -Dq 2 (l + rDq 2 + ■ ■ ■) (63) 
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which is definitely negative if r > 0. This is the usual situation where the 
waves falls off exponentially in time. If r < then the first two terms of the 
expansion cannot be trusted for q 2 > 1/\tD\. But the first two terms of the 
expansion cannot be trusted for larger values of q even if r was positive. 

Now consider real, positive uj and complex q. The shear dispersion rela- 
tion is 

q 2 = ^( i + T LU + ---) . (64) 

The wavenumber can be decomposed into its real and imaginary parts as 
Q = Qr + iQi- The dispersion relation to the required order is determined by 
the solution to the following equation: 

q R -Q 2 i-^- + i[2<lR<li-p) =0. (65) 
Obviously q^qj = u/2D > 0. The shear wave behaves as 

$ = ijj e i(q R z-"t) e - qi z _ (gg) 

For a wave moving in the positive ^-direction, qn > and qj > 0, and 
the wave is damped in space. For a wave moving in the negative ^-direction, 
qn < and qj < 0, and the wave is again damped in space. Hence, these shear 
waves are damped in space for any finite value of the relaxation — oo < r < oo 
no matter what its sign. 

We are not aware of anything unphysical about a negative shear relax- 
ation time according to the definition of it given in this paper. See also the 
analysis and discussion of this point in [25] . The only unpleasantness associ- 
ated with a negative relaxation time may be numerical instabilities that could 
arise in relativistic viscous fluid dynamics due to finite numerical resolution. 

7 Conclusions 

In this paper we have used the gauge/gravity correspondence at finite tem- 
perature to study the properties of the shear mode relevant to a four dimen- 
sional world. We confirmed a general formula for the shear diffusion constant 
- equivalently the dimensionless ratio of shear viscosity to entropy density. 
This quantity may be computed if the five dimensional metric is given. It 
reproduces all known results in the literature. The method we use allows us 
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to determine the shear dispersion relation to arbitrary high order in the fre- 
quency and wavenumber. In particular, the main result of this paper was to 
obtain a general formula for the coefficient of the next order term in the shear 
dispersion relation. It has the dimension of time, and it is often referred to 
as a shear relaxation time, although that is somewhat of a misnomer as we 
explained earlier. It too reproduces all known results in the literature. We 
applied these formulas to models of confinement in QCD, namely, the hard 
and soft wall models. These results may be useful in relativistic viscous fluid 
descriptions of high energy nuclear collisions at RHIC and LHC. 

We are currently pursuing several avenues of research. We would like to 
perform the same type of analysis for the sound mode. However, that is a 
much more difficult undertaking, and seems to require more information than 
is necessary for the shear mode. A precise criterion for determining which 
theories these results apply to has not yet been given in the literature. We 
are investigating the coupling of arbitrary fields to the gravitational field in 
an AdS space to address both of these problems. 
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Appendix 



The Riemann-Christoffel tensor Rap^ for the background metric has the 
following independent components. All the rest are either related to these 
by the algebraic properties of Rap^u or are zero. 
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(67) 



The diagonal elements of the Ricci tensor R^ v = g a/3 R a ^/3 U are nonzero 
while the off-diagonal ones are zero. 
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The curvature R = R x x is given as follows: 

R — Fq + 3F X + F r . 
Lastly we come to the invariant 

RjxvXcyR^ = g^g 13 " g Xl g aS Rap-fsR/ivXo 

For the background metric this reduces to 
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which can be written as 
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Figure 1: The ratio of the shear viscosity to entropy density as a function of 
temperature in the soft wall model. 
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Figure 2: The shear relaxation time as a function of temperature in the soft 
wall model. The dashed curve is the high temperature expansion given in 
the text. 
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